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In this paper we study the time evolution of a class of two-level systems driven by periodic fields 
in terms of new convergent perturbative expansions for the associated propagator U(t). The main 
virtue of these expansions is that they do not contain secular terms, leading to a very convenient 
method for quantitatively studying the long-time behaviour of that systems. We present a complete 
description of an algorithm to numerically compute the perturbative expansions. In particular, we 
applied the algorithm to study the case of an ac-dc field (monochromatic interaction), exploring 
various situations and showing results on (time-dependent) observable quantities, like transition 
probabilities. For a simple ac field, we analised particular situations where an approximate effect 
of dynamical localisation is exhibited by the driven system. The accuracy of our calculations was 
tested measuring the unitarity of the propagator U(t), resulting in very small deviations, even for 
very long times compared to the cycle of the driving field. 
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I. INTRODUCTION 

Periodically (or more generically quasi-periodically) 
driven quantum two-level systems are of basic impor- 
tance in many physical applications an have been widely 
studied since the pioneering works of Rabi Q , of Bloch 
and Siegert Q and of Autler and Townes || (see also ||] 
for more recent discussions and for other general 

references on the subject). The time evolution of such 
systems is governed by the Schrodinger equation (we 
henceforth adopt h = 1) 



i- |*) = H(t) |*) 



(1) 



where |*) = = and H(t) is the Hamilto- 

nian of the system. We may consider, for instance, a spin- 
1/2 system under the influence of a time-dependent (pe- 
riodic) magnetic field B(t). In this situation, H(t) takes 
the usual form H(t) = — ^B(t) ■ cr, where a — (<ri, o~2, 03) 
are the Pauli matrices. The interest in the solution of (|l|) 
in this case is not restricted to the investigation of the 
quantum system. As first pointed by Feynman, Vernon 
and Hellwarth jio) (see also the discussion in [jy]), the 
quantum system is equivalent to the classical Hamilto- 
nian system describing a classical gyromagnet precessing 
in a magnetic field: j^S = —B{t) x S, where S is a (three 
dimensional) unit vector. In fact, the methods described 
below can be directly applied to the analysis of this sys- 
tem as well, since the components of S can be written in 
terms of the components ip±(t) of |^) M,0. 

The evolution of the systems governed by (jl|) is de- 
termined by the time evolution operator U (t, s) (also 
known as the propagator) which connects the state |^(s)) 
at time s to the state | at time t: \^(t)) — 
U(t, s) |*(s)) . Defining U(t) = U(t,0) one has U(t, s) = 



U(t)U(sy . For a time-dependent Hamiltonian the prop- 
agator U(t) can be computed via the Dyson expansion: 



H{ti)dt x 



tn-l 



H{t n ) dt r , 



(2) 



Although (|) gives a straightforward manner to compute 
U(t), the series in the r.h.s. is not generally uniformly 
convergent in time. For practical purposes this gives rise 
to difficulties when one is interested in the large-time be- 
haviour of the system. For instance, if one considers a 
periodic Hamiltonian of the form H{t) = J2 m H m e muJt , 
with Ho ^ 0, two successive integrations in (||) would 
produce a linear term in t. Higher order terms in t 
would appear with further integrations. These polyno- 
mial terms are known as secular terms and they plague 
the expansion of U (t) in such a way that its uniform con- 
vergence is spoiled. 

Of particular interest is the situation where the 
Schrodinger equation (|l|) takes the form 

i^|*> =fli(*)|¥) , with H x {t) := eaa - , 

(3) 

where f(t) is a periodic function of time t and e is a real 
constant. By a rotation of 7r/2 around the 2-axis, we get 
the equivalent system 



ij t |*> = H 2 (t) I*) 



with H 2 {i) := eo-x + f(t)a 3 , 



(4) 



where |$) = |$(t)) = exp(-i7rcr 2 /4) | *(£)). The Hamil- 
tonian in (^|) can be interpreted as describing a system 
with an unperturbed diagonal Hamiltonian Hq := e<73, 
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subjected to a (periodic) time-dependent perturbation 
Hi(t) := —f(t)ai, inducing a time-depending transition 
between the unperturbed eigenstates of H . Of course, 
(||) can be also interpreted as a spin- 1/2 system under a 
magnetic field B(t) = (2/(i),0, -2e). 

In we investigated the system described by (||) or 
(Q) in the situation where / is a periodic or quasi-periodic 
function of time and e is "small" and a special perturba- 
tive expansion (power series in e) was developed, whose 
main virtue is to be free of secular terms. The algo- 
rithm employed involves an inductive "renormalisation" 
of a sort of effective field introduced through an expo- 
nential Ansatz (the function g, to be introduced below). 
Moreover, in the periodic case, we have been able in |7]]|] 
to prove rigorously the uniform convergence of our ex- 
pansions for all t, for e sufficiently small, provided the 
function / satisfies the special conditions I, II and III we 
describe below. 

A feature of our method is the fact that we are able 
to present the complete (and convergent) e-expansion for 
the secular frequency ft (also known as Rabi frequency) 
associated to the solution of (H|)-(||). This is particularly 
important for the qualitative study of the large-time be- 
haviour of that solutions. Another relevant conclusion is 
that, under conditions I-III, our method provides a com- 
plete perturbative solution of (||)-([I|) for the monochro- 
matic field (ad-dc field), except perhaps for spurious sit- 
uations Q. 

Since our series are uniformly convergent in time, it is 
possible to use them to study the long-time behaviour of 
observable quantities in a controlled way. It turns out 
that our expansions are also very practical for numerical 
calculations, providing very small errors even with few 
terms. 

In this paper we describe the algorithm employed in 
the numerical computation of our perturbative solution 
and show the results obtained in the particular case of 
the monochromatic interactions. We begin in Section |fi| 
by given a quick review of the perturbative method de- 
veloped in Section [II presents a description of 
the numerical algorithm developed to compute the ex- 
pansions. Next, in Section |v|, we show and discuss the 
results obtained in the case of monochromatic interac- 
tions. Finally, in Section [v| we draw some conclusions 
and final remarks. 



II. DESCRIPTION OF THE METHOD 

In B , it was shown that the propagator U(t) associated 
with the system (||) can be written as 



U{t) 



R(t)(l + ig S(t)) -ieR(t)S(t) 



-ieR(t)S(t) R(t)(l - mS{t)) 



, (5) 



where 



and 



R(t) := exp 



S(t) 



•* / (f(r)+g(T))dT 



R(t)~ 2 dr , 



(6) 



(7) 



where g is a particular solution of the generalised Riccati 
equation 



g(t)' - ig{t) 2 - 2if(t)g(t) 
and #o = 5(0)- To solve (^|) we pose 

oo 

g(t) = J2G (n) (t)e n , 

n=l 

where G^(t) := q(t)c n (t), with 



it 2 = 0. 



(8) 



(9) 



q(t) := exp J f{r)dT 

Inserting (|) into (|) yields a sequence of recursive equa- 
tions for the coefficients c„, whose solutions are 



ci(i) = a\q(t) , 
ca(«) = q(t) 

c n {t) = q{t) 



i I (a 1 g(r) - q(r) ) dr + a 2 
o 

i I / y~] c p {T)c n - p (T) dr + a 
V° p=i / 



(10) 
(11) 

(12) 



for n > 3, where the a n are arbitrary integration con- 
stants. Our strategy consists in fixing these constants 
inductively to cancel the secular terms. For instance, to 



cancel the secular term in C2, the integrand in (11) can- 
not contain a constant term. Defining the mean value of 
a quasi-periodic function h by 

M(h) := lim -L f h(t)dt, 



we must require that M(a\q 2 



0. This leads to 



a\ = M(q 2 )/M(q 2 ), provided M(q 2 ) ^ 0. We remark 
that the mean value of a quasi-periodic function h equals 
its zero order Fourier coefficient. In [^]-|| we identified 
three classes of quasi-periodic functions / for which this 
procedure of elimination of secular terms can be applied 
to all orders (i.e., to all functions c„). These classes are 
defined by the following conditions: 

I. M(Q ) ^ , 

II. M(Q ) = but Af(Qi) ^ 0, 
III. M(Q ) = M(Gi) = but M(Q 3 )^0, 
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where Q (t) :=q(t) 2 , 

Qi(t) := Qo(t) / (Qoiry 1 - M (Qq l )) dr , (13) 
Jo 

Qs(t) := Qo(t) f (Qi(r) - M(Qi)) dr . (14) 
Jo 

Under condition II, the Ansatz (0) has to be slightly 
modified to 



fl(t) =J]5< n )(t)A B , with A:=e 2 

n=l 

and G"(i) := q{t)e n (t). The solutions for e„ are 



d(t) =g(t) ( -i J q{T)- 2 d T + 0, 



(15) 



Sn(*)=9(0 



(16) 
(17) 



for n > 2. Under condition III, we found in || the same 
solution for g as in I, namely (S). The difference between 
the solutions obtained in I and III lie in the set of con- 
stants a n , which is conventionally fixed to remove the 
secular terms in each case. 

In pf] we have shown that the e-expansion (|j]) or jl5| ) 
converges uniformly for |e| small and for all t, provided 
/ is periodic and satisfies I or II. In H we extended this 
result for case III. As we will discuss soon, case III is par- 
ticularly important for the investigation of the dynamical 
localisation effect (also called, less properly, coherent de- 
struction of tunnelling) for monochromatic interactions. 

The secular frequency fl, in conditions I or III, is given 

by §§ 



Q = M(f) + M(g) = F + ]T M e n 

n=i 

and, in condition II, by 

oo 

V = Fo + ^M (& n) ) A". 



(18) 



Hence, with our previous definitions, for condition I, we 
have 



n = F + eaiM(Qo) 

+ e 2 [ialM(Q 2 ) - tM(Ci) + a 2 M(Q )] 
+ e 3 [2a x M (fi 3 ) + a 3 M{Q )} + 0(e 4 ) , 



(19) 



where 



e 2 (*):=Qo(*) / (Q (T)-M(Qo))dT. 
Jo 



One easily sees, by computing the zero order Fourier co- 
efficient of Q 2 (details in ||), that M(Q 2 ) =0 whenever 
M(Qq) = 0. Hence, for condition II, we get 



(20) 



n = F - ie 2 M(Q 1 ) + 0(e 4 ) 
and, for condition III, 

O = F + 2a 1 e 3 M(Q 3 ) + £>(e 4 ) . (21) 



III. NUMERICAL COMPUTATION OF THE 
SOLUTION 

We present here a step-by-step algorithm which en- 
ables us to compute the propagator (^) solely from the 
Fourier coefficients of the (periodic) interaction /. This 
algorithm is based on the convergent expansions pre- 
sented in 0, which expresses the propagator U(t) in 
terms of its Floquet form. 

Let us first begin with some conventions. We sup- 
pose that the interaction function / is periodic with 
period T w = 2-k/oj and that its Fourier decomposition 
/ = J2 n F n e muJt has only a finite number of terms. 
Since / is real, excluding the constant term Fq, f must 
have an even number of non-vanishing Fourier coeffi- 
cients, say 2 J, with J > 1. Denoting the set of integers 
{neZ,n^0\F n ^0} by {n%, . . . , n 2 j}, we may write 

2J 

f(t) = F + J2fae in « ut , 

a=l 

with the convention that n a — — n^j— +i> f° r all 1 < a < 
J, and with f a = F„ a . Clearly f a = f 2 j- a +i, I < a < J. 



A. Interactions with Fq — 

In this case, the Fourier decomposition of the functions 
q and q 2 can be written as 



and q(t) 2 = ^ Q 



(2) ^imuit 



mG2 



(2) 

The coefficients Q m and Q m are of basic importance in 
the numerical computation of the propagator (^J). They 
can be obtained in a closed form from the coefficients F m 
of /. Explicitly, we have 



2.1 



— P l 1f 



E <K p >™)II 



Pi, ■■■ ,P2J=0 



a=l 



J_ (A. 

p a \ \n a uj 



(22) 



for all to, where 



2,1 



P = P(Pl,--- ,P2j,n!,... ,n 2J ) :=^Piiii,, 

6=1 
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and 



If 



2J , 



The symbol 6(P, m) denotes the Kronecker delta: 
8(P,m) = 1, if P = m and 6(P,m) = 0, otherwise. To 
compute the Fourier coefficients of q 2 , we simply note 
that q 2 is obtained from q by the substitution / — » 2/. 
Hence, for all m, 



2.7 



g 



(2) 



Pi; 



■■ ,P2J=0 



0=1 



1 



2/g 



(23) 



Formulas ( p2| ) and ( p3|) can be computed either analyti- 
cally (in some cases) or numerically. For the monochro- 
matic interactions, a closed form in terms of Bessel func- 
tions of first kind is obtained. The exact result will be 



presented in Section IV. 

(2) 

Once wc know the coefficients Q m and Q m , we pro- 
ceed to compute the particular solution (||) or ( [l5| ) of 
the generalised Riccati equation (||). To decide whether 
we use @ or Jig), we must check which condition I, 
II or III is satisfied by /. Hence, we have to look at 
the mean values of Qq, Q\ and Q3. One obviously has 
M(Q ) = M(g 2 ) - Q { o } . To obtain Af(Qi) and M(Q 3 ), 
wc need the Fourier decompositions of Qi and Q3. From 
the definitions ( fl3| ) and (14), after some simple compu- 
tations, we get 



M(Qi) = - 



Q { -l(Q 



(2) 




,(2) 



m ^0 



7/1 



and 



Af(Q 3 



E 



)(2) 
> — m 



cj^ ■« — ' nm 



(2) 



V0 ^n-m 1 



The numerical value of M(Qi) and M(Q 3 ) can be calcu- 
lated trivially from the above expressions once Q m are 
known. 



1. Computing g(t) in Cases I and III 

We remember that condition I applies whenever / sat- 
isfies M(Qo) 0. In this case, by properly fixing the 
constants a n , we can completely eliminate the secular 
terms from the functions c„ For this reason, wc 

may write 



Cn(t) = E 



The Fourier coefficients Cm are obtained from equations 
(|io|)-(|l2|). Their inductive structure is given by the rela- 
tions 




a2<3m , 



n-1 



E^" 



n-p) 
"2 



p=l 



+ CKnQm , for Tl > 3 . 



The constants av 



ficients and Cm' , for p < n 



have closed forms in terms of the coef- 
1. Since they involve 
somewhat large expressions, we refrain from writing them 
here. The complete expressions can be found in ||. We 
may see that the whole inductive structure of the coeffi- 
cients Cm^ is known and, therefore, the computation of 
c n is just a matter of numerically evaluating the above 
expressions. 

Since in case I g is given by (^|), we may write 



,(t) = q(t) E c n (t)e n =: E E G^e w e r 



n=l 



n=l \m£2 



(24) 



where Gm' is given by the convolution 

m — / y ^im-p^p j 

pSZ 

whose numerical value can be easily computed since we 

(n) 

already know Q m and C m , for all m and n. This 
gives g for all i under condition I. For condition III 
(M(Q ) = M(Qi) = 0, but M(Q 3 ) ^ 0) one has es- 
sentially the same solution except for the constants a n 
which are calculated differently from condition I. Their 
formulas can be found in [pi. 



2. Computing g(i) in Case II 

Condition II applies whenever / satisfies M(Qq) = 0, 
but M(Qi) ^ 0. In this case, the perturbative solu- 
tion of the generalised Riccati equation (||) is given by 
(|l5|). The constants (3 n which appear in the functions e„ 
(equations (|l|)-(|l7])) can be chosen such that no secular 
terms emerge ||. Hence, with this particular choice of 
the constants /?„, we may write 



e n (t) = E E^e* 



The recursive structure of the coefficients E m ^ is given 

by [§ 
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(1) 



E 



E Qm-nQ-n p. 
mo 

y 



A + E 



Qrn — ni—r 



E^M 



n-p) 



-<?« E E**M 



n-p) 



1 



(ni + 712)0; 



+ /? n Qm 7 for n > 2 . 

The constants /3 n have also closed forms [|| in terms of 
the known Qm and Em , for all to and p < n — 1 . Thus, 



the solution ( |15|) can be evaluated numerically in any 
order for all t. For future convenience, we write 



1® = q (t) J2 e»(*)A n = : E E 1 A " 



n— 1 \m€2 



where Gm' is given by 

pez 

which completely specifies g for all i under condition II. 

3. Computing the Propagator in Cases I, II and III 

Once the coefficients Gm' (for cases I and III) or G$ 
(for case II) are known, the propagator U(t) expressed 
in (||) can be computed in a straightforward manner, as 
will be shown now. We illustrate our procedure with the 
coefficients Gm of conditions I and III. For condition II, 
where the coefficients of g are Gm , the discussed proce- 
dure has to be adapted with self-evident modifications. 

We begin by defining 



G m (e) 



(26) 



The secular frequency ( |18[ ) is clearly given by (see (24)) 

oo 
n=l 

since we are supposing that F = 0. Hence, by (|24|), 

g(t) = n+J2G m (e)e im " t . (27) 

Looking at expression (g) for the propagator, we see that 
the Fourier series of R (see (0)) can be computed if we 



first find the Fourier decomposition of 



W(t) := exp ( -i j g(r) dr 



Indeed, since R — qW, we obtain the Fourier series of R 
by taking a convolution of the coefficients of q and W. 
It is easy to see that 



W(t) 



expf-^TiW™' j 

\ mGZ / 



(28) 



with 



H m = H m {e) := , for to ^ 

TOCJ 



and i?o = 0- Moreover, 



(25) Writing 



7/0) :=* E H " 



mez 



(29) 



and using (|28j), we find 



x flni ' ' ' H np H m - Np , for to ^ 



and 



/ OO ,_ -|Np+l 

WbsWb(e)=e-*r/W l + J^ 



ri\ , . . . , n p GZ 



with 7Vp := X)a=i n a- Using now ( |29| ) and the fact that 
i? = qW, we conclude that i? can be written as 



R[t) =e Mm e 



(30) 



with the coefficients R m given by the convolution 

Rm = E *5 

p—rn 

'W p . (31) 

p€Z 

This finishes with the computation of R in terms of its 
Fourier series (p0|). We note from formulas (|2^)-(|3^) that 
once the coefficients Gm' of g are given, we can numeri- 
cally evaluate G m , H m , W m and, hence, R m with a trivial 
computer code. 
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Next we proceed to compute the Fourier series of S 
(see First we find the Fourier coefficients of Rr 2 . 

This is an easy task since R~ 2 is obtained from R by 
replacing (/ + g) — > — 2(/ + g). Hence, we must replace 
H m —2H m and q — > q 2 . Consequently, we get 



R(t)- 2 = e 2int R(-V e m ^. 



(32) 



where 



r(-v = Vo (2) w^ 2 

m / ^ ^ m—p r T p 



to Q m and Q„ , derived from the Fourier coefficients of 
the interaction function /. In short, to help us visualise 
the necessary steps towards the computation of U(t), we 
may draw the following "chain" : 



/(*) 



w m ,wL 2) 



n 



) o (2) 



r<( n ) „ 
, go 



a 



(n) 



R{t),S(t) -» U(t). 



with 



2 p+i 



W m ^ ee W^(e) = e 2 ^ ^2H m + £ (p , 



2J H ni ---H np H, 



m-N„ 



, for m ^ 



ni , . . . , n p £2 



and 



IF, 



(-2) 



2 p+i 




Til,... 1 1p£ 



Now 5 is obtained by a simple integration of R 2 . A 
trivial computation from (p2), gives 



S(t) =( 7 + e 2 ^^S m e 4 ' 



(33) 



B. Interactions with Fo 7^ 

When F 7^ 0, one automatically has M(Q ) = 0, ex- 
cept perhaps when 2Fo = kuj, for some integer k. These 
facts were shown in Since the situation where 

2Fq = kuj, for some integer fc, was nowhere investigated 
in our previous works we will ignore this possibil- 

ity by henceforth assuming that 2Fo 7^ ku>, for all k £ Z. 
Since M(Qq) = 0, condition I is never satisfied when 
F) 7^ 0. The convergence of our expansions for Fo 7^ in 
condition III, however, has not yet been studied. Hence, 
we will only consider condition II. In this case {tJ, the 
function q turns to be 



q(t) = e tFot Q m e imui 



and the functions e„ (see (|l6|)-(|l7|)), 



e n (t) = e 



-iFat 



(36) 



(37) 



with 



R 



(-2) 



muj + 2ft 



and ctq := — S n 



(34) 



III fzl 



We assume that muj + 2fl 7^ for all m S Z (see the 
discussion of crossings in (7)). 

We have found expressions for F and 5 in terms of its 
Fourier series. This series converge absolutely and uni- 
formly as we showed in |t]|| . To compute the propagator 
given in (||), we still need go, which can be easily obtained 
from (127 



g a = ,g(0) = fl + G ™( e ) 



(35) 



,,.#0 



Formulas (pOj) , ( |33[ ) and ( |35| ) can now be used to eval- 
uate U(t) for all times. We stress that these formulas 
depend essentially on the Fourier coefficients Gm^ of g. 

(n) 

These, in turn, depend on C m which are direct linked 



When F 0, the coefficients in the above expres- 
sion assume a simple form J7J: 



F« 



EQm-\-aQa 
CLLO + 2F ' 

n-1 ^ 



Etf = E E 

p=l a,b£Z 



) t p(p)p(™-P) 

(a + 6)w - 2F 



for n > 2 . 



Finally, from (|36|) and (|37j), we conclude that g given in 
( p5| ) assumes the form 

where 



G m = G m (e) 



E g ™ )a "' 

n=l 



G 



with A 



and 



This gives g in terms of its Fourier series. 

To compute the propagator U(t) accordantly to (||), 
all we have to do is follow the procedure detailed in Sec- 



tion 

e — » 



[II A 3 with the replacement of G 



(n) 



Gm^ and of 



A. 



IV. NUMERICAL RESULTS FOR THE 
MONOCHROMATIC INTERACTIONS 



We now apply the algorithm discussed in Section f 
to study the system ([|)-(|]) under the influence of 
monochromatic interactions (ac-dc field): 



f(t) = F + ipcos(ujt) 



(38) 



which are relevance for many physical applications. With 
the conventions introduced in Section III , we have J = 1 , 
fi = h = <f/2, ni : 
formula (|23|) gives 



-ri2 = — 1. A simple application of 



Qo(t) 



riGZ 



where J n is the Bessel function of first kind and order 
n and where we defined xi := and X2 := 2Fo/u. 

Depending on the parameters xi an( l X2j the function / 
given in ( |38| ) satisfies one of conditions I, II or III, except, 
perhaps, for spurious situations. The detailed analysis of 
these facts can be found in Ji| . Table | summarises the 
conclusions presented in || and gives a classification of 
the conditions satisfied by / as a function of the param- 
eters Xi and Xi- 



Label 


Xi 


X2 


Condition 


(A) 


not a zero of J m 


— m £ Z 


I 


(B) 


zero of J m 


-m e z 


III 


(C) 


any 


not integer 


II 



TABLE I. Classification of the conditions satisfied by the 
monochromatic interactions f(t) = Fo + <p cos(o;t) as a func- 
tion of the parameters xi = 2<p/uj and \2 = 2-Fo/w. We 
labelled by (A), (B) and (C) the three possible cases. 



Although not indicated in Table [| there is in case (C), 
on each interval (k, k + 1), k = 1, 2, . . . , a special value X2 
of X2 (depending on xi) f° r which M(Qi) = 0, and we 
would be out of condition II |8) . We refrain from study- 
ing this rather spurious situation here. A more detailed 
analysis of this case can be found in M . 



We next show some graphical results of transition 
probabilities calculated via the algorithm described in 
section pi] for the situations (A), (B) and (C) described 
in Table|[ Let |$ + ) = (J) and |$_) = (°) be two or- 
thogonal states of a system described by ([I]) (the eigen- 
states of the unperturbed Hamiltonian H of Hi in (ph). 
The probability for the transition from the initial state 
|$+) to the final state |$-) at time t is given by 

P(t):=\{$ + \U(t)\$-)\ 2 = \U 12 (t)\ 2 . (39) 

We can evaluate P(t) numerically using the methods de- 



scribed in Section III to compute U(t). To estimate the 
accuracy of our calculations, we tested the unitarity of 
the time evolution operator, U(tyU(t) — 1, and consid- 
ered the quantity 



N(t) :=|tfii(t)| a + |EM«)| S 



1 



which should be identically equal to for unitary U(t). 

Let us first consider case (A) of Table | with u = 1.0, 
Xi = 2/cj and xi = Fq = 0. Figure |l| shows graphs of 
P{t) and N(t) for e = 0.01, e = 0.10 and e = 0.40, plot- 
ted from t = to t = Tq = 2-k/VL : in units of T u = 2tt/uj, 
the cycle of the external field. The calculations were 
performed using an expansion for g up to 0(e 25 ). We 
took all the Fourier coefficients (generically called F m ) 
involved in the computations of U (t) within the range 
to = —40, . . . ,40. From the deviations of N(t) from 0, 
we can infer very small errors in the calculations, lead- 
ing to very accurate values of P(t). For e = 0.01 and 
e = 0.10, we have errors of the order of only 4.0 x 10~ 5 %. 
For e = 0.40, the errors jump to 6.0 x 10~ 2 %, indicat- 
ing that the parameter e is coming close to the radius 
of convergence of our expansions. Since in case (A) we 
are under condition I, the secular frequency $1 is given 
by ©, hence il = 0(e) (we choose F = X2_ = 0). For 
e = 0.01, from the full expansion of Q (see (|18|)), we get 
T n = 450T W . For e = 0.10, T n ^ 45T W and for e = 0.40, 
Tq = 25T W . We may notice from the graphs of Fig- 
ure 0, that the transition probability P(t) behaves like a 
Rabi oscillation with a frequency f2, i.e. P(t) = sin 2 (fit). 
Small oscillations of frequency lu, with amplitudes of or- 
der |e|, are superposed with this Rabi oscillation, leading 
to a quasi-periodic evolution for the system (as ensured 
by the Floquet theorem). 

We now consider case (B) of Table Q with u = 10.0, 
Xi = Xi, x± being the first positive zero of Jo, and 
X2 = Fq = 0. Since we are under condition III and 
F = the secular frequency becomes, according to (pi]), 
£1 = C(e 3 ), a fact first pointed in ||. This weak depen- 
dence on e indicates long transition times for the proba- 
bility amplitude P(t). This phenomenon is known as the 
dynamical localisation effect [[l2| (sec || for a general 
criterion and for a more complete list of references). To 
compute the solution in this case, g was expanded up to 
e e . This contrasts with the expansion for g in the previ- 
ous situation, where g has been computed up to C(e 25 ). 
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The reason for this lies in the fact that under condition 
III the expressions for the constants a n (see (|ic|)-(|l^)) 
are somewhat more intricate than those of condition I. 
Fortunately, however, the expressions for the constants 
a>i, a 2 and are not so complicated and this suffices 
to give an expansion of g free of secular terms up to the 
e 6 term (see B for details). If we limit ourself to study 
situations of small e, i.e. well inside of the radius of 
convergence of our expansions, excellent results can be 
achieved. To improve the accuracy of our results, it is 
better to use a larger value of w. 

Figure shows graphs of the numerical values of P(t) 
and N(t) for e = 0.01, e = 0.10 and e = 0.20, plotted 
from t = to t = Tq = 2tt/CI, in units of T u = 2tt/oj. 
The errors in Figure ^| (measured out of the deviations 
of N(t) from zero) are greater than those presented in 
Figure |l| due to the lower order e-expansion of g in the 
former case. However, for small e, the errors are still 
small, being of the order of only 3.0 x 10 -3 % for e = 0.01 
and of the order of 3.0 x 10 _1 % for e = 0.10. For larger 
e (= 0.20), the errors reach the value of about 1.0 %, 
indicating that a higher order e-expansion is needed to 
improve accuracy. We observe that the plots of P(t) in 
Figure @ have the same qualitative aspect of those in 
Figure nl There is a predominant Rabi oscillation (with 
frequency £1) aspect superposed by minor oscillations of 
frequency tu whose amplitudes are possible bounded by 
|e| 3 terms. This fact, however, has yet not been proven 
directly from our expansions, so it has to be faced more as 
a qualitative analysis rather than a rigorous quantitative 
one. 

The main distinction between the graphs of Figure |l| 
and those of Figure |2| is, undoubtedly, the Rabi oscillation 
period Tq. In the later, the long time needed for the sys- 
tem to transit from the initial state |$+) to the final state 
1$-), compared to the basic cycle T u of the external per- 
turbation, is the (approximate) dynamical localisation 
effect. The secular period Tq obtained for the situations 
studied in Figure | were Tq = 1.6 x 10 9 T u for e = 0.01, 
Tq 1.6 x 10 6 for e = 0.10 and Tq ^ 2.2 x 10 5 T w for 
e = 0.20. These are somewhat much larger values of Tq 
than those of Figure |l|. 

As we have mentioned, the graphs of Figure || were 
computed with \i — X U where x\ is the first positive 
zero of Jo- It is interesting to test our solution consider- 
ing other possible zeros of Jo. In Figure || we show plots 
of P(t) and N(t) calculated with \2 — the second 
positive zero of Jo (we also took uj = 10.0 and, of course, 
\2 = 0). The three situations presented in Figure ^ cor- 
respond to e = 0.10, e = 0.20 and e = 0.30. We may 
note that the qualitative behaviour of P(t) presented in 
Figure |^ is the same of Figure |^. In particular, the ef- 
fect of dynamical localisation is preserved since we still 
have fl = C(e 3 ) when xi — x i- We may note from the 
deviation of unitarity N(t) shown in Figure || that the 
errors, compared with the ones in Figure | for the same 



values of e, are smaller by a factor ~ 1/10. We may un- 
derstand this fact as follows: since X2 > x±, the strength 
ip of the interaction presented in Figure || is greater than 
in Figure ^, hence, for the same values of e, the effective 
perturbation e/ip is smaller in the former case, leading to 
a more precise perturbative computation. 

We now investigate case (C) of Table Bj. We consider 
u> = 1.0, xi = 1.0 and \2 =0.3. This implies condition II 
and we have to follow the prescriptions of Section III B| . 
Figure ^ shows the results obtained for P(t) and N (t) for 
three values of e: 0.05, 0.10 and 0.20. The various plots 
were calculated using an expansion of g up to e 20 . It was 
not necessary, thus, to use large u> to ensure convergence 
of the expansions, as we did in the situations presented 
in Figure § (condition III). As usual, P(t) and N(t) were 
computed from t = to t = Tq in units of T u . The 
qualitative behaviour of the transition probability is sig- 
nificantly different from the previous results (Figures |l| 
and |). This is a consequence of the non-vanishing con- 
stant field Fq presented in the interaction /. Indeed, the 
secular frequency f2 is now given by ( po| ) and, hence, is 
of order of Fq. For e = 0.05, we obtained Tq = 6.3T W , for 
e = 0.15, T n = 6.0T„ and for e = 0.20, T n 5.0T W . Since 
we choose Fq as the same order of to, there is a strong 
competition between the Rabi oscillation (governed by 
fi) and the external field oscillations. This leads to the 
patterns shown in Figure ||, which do not behave purely 
like sin 2 (fW). In particular, P{t) < 1 for all times, lead- 
ing to the conclusion that the state |$+) never transits 
completely to |$_). One sees, moreover, that the tran- 
sition amplitude P(t) tends to zero as e — * (c.p. fl4l|), 
below), much in contrast to the cases pictured in Figs. 
[I], |^ and |^. The leading e-dependence of P(t) can be al- 
gebraically determined, in principle, but this was not yet 
performed due to the complexity of our expansions. 

It is interesting to note from the graphs of Figure [| that 
the transition probability gets closer to 1 as e increases. 
It would seem that for e large enough, we could have 
P(t) > 1. However, our expansions would not converge 
in this case, since e would be greater than the radius of 
convergence. To illustrate this situation, let us consider 
the trivial case where Xi = and \2 — 2Fq/lu is not an 
integer. Since now f(t) = Fq = LUX2/2, the Schrodinger 
equation ([|) becomes time- independent, with a simple 
Hamiltonian given by 



H 2 = ecri + F cr 3 



F 
c 



c 

-F 



The propagator U(t) can be computed by elementary 
methods (f.L, by diagonalising the Hamiltonian), lead- 
ing to 



U(t) = cos(uj t)i - ( iFQ(X3 + e(Ti - 

UJ 



(40) 
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where lu :— y/F^ + e^. Hence, 
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P{t) 



F 2 



■ sin 2 (uJot) . 



(41) 



Note that P(t) < 1 for all e and all F . For f(t) = F , 
the generalised Riccati equation (^) admits a particular 
solution given by the constant go = —Fq + \J F 2 + e 2 . It 
was shown in Q that the e-expansion ( |i~5| ) for f(t) = F a 
coincides, as expected, with the Taylor expansion (cen- 
tred at e = 0) of go. Thus, our method has a clearly 
restricted region of convergence defined by |e| < \F \. 

Let us see what these last considerations mean numeri- 
cally. First we set Fo — 0.4. Figure [| shows plots of P(t) 
and N(t) for three critical values of e: e — 0.30 < Fo, 
e = 0.40 = Fo and e = 0.43 > Fq, where our expan- 
sions are not supposed to converge. We also show in 
Figure pi (dashed lines) plots of P(t) calculated accord- 
ing to (}4l|). Looking at N(t) and at the deviation of the 
perturbatively computed transition probability from the 
one calculated via ([4l|), we conclude that the region of 
convergence of our expansions is restricted to |e| < |Fo|, 
as expected. It is important to note that when e < Fq, 
our computation of the propagator matches exactly with 
the non-perturbative solution (fiof). This can be seen 
from comparing P{t) calculated perturbatively and non- 
perturbatively via ([4~ij), as shown in Figure |5|. 
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We stress that, at least for the case of monochromatic 
interactions examined in Section IV, the errors obtained 
are very small and bounded as time increases (see the 
behaviour of N(t) in figures of Section IV). This is due 
to the absence of secular terms in our perturbative ex- 
pansions and its uniform convergence in time. As a con- 
sequence, one can study the long-time behaviour of the 
quantum system (||)-(|]) in a controlled way. 

Another important feature of our method is that it 
can be numerically implemented with relatively simple 
computer codes. Indeed, the algorithm described in Sec- 
tion III to calculate the unitary propagator U (t) consists, 



essentially, in simple computations involving the Fourier 
coefficients of the functions q and q 2 , which are known in 
a closed form (see (|||) and (|3|)). This is an important 
advantage against other perturbative approaches, based 
f.i. on the Dyson expansion (|^) which, in general, cannot 
be evaluated in a simple manner (not to mention the fact 
that such expansion is not uniform convergent in time, 
as we have stressed in the introduction of this paper) . To 
sum up, our method is not only mathematically rigorous, 
but also very useful for practical purposes, where it can 
be applied with great generality, leading to very accurate 
results. 
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FIG. 1. Plots of the transition probability P(t) (left col- 
umn) and the deviation from unitarity N(t) (right column) 
as a function of time (measured in units of T u ) for various e. 
We considered case (A) of Table |, with lj = 1.0, xi = 2/w 
and X2 = 0. We used e = 0.01,0.10,0.40 in the top, middle 
and bottom rows, respectively. 




FIG. 2. The same of Figure |l|, but we now considered 
case (B) of Table Q, with lo = 10.0, xi = ^l, x i being the first 
positive zero of Jo and \2 = 0. We used e = 0.01, 0.10, 0.20 in 
the top, middle and bottom rows, respectively. Note that the 
time scale is multiplied by the factor in the square bracket. 




FIG. 3. The same of Figure g, but now with xi = x 2, 
%2 being the second positive zero of Jo and \i = 0. We 
used e = 0.10, 0.20, 0.30 in the top, middle and bottom rows, 
respectively. 




1.0 2.0 3.0 4.0 5.0 1.0 2.0 3.0 4.0 5.0 



FIG. 4. The same of Figure hi but we now considered 
case (C) of Table |, with u = 1.0, xi = 1-0 and X2 = 0.3 (not 
integer). We used e = 0.05, 0.10, 0.20 in the top, middle and 
bottom rows, respectively. 
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FIG. 5. The same of Figure but now with \i = and 
X2 = 0.8 (Fo = 0.4). We used e = 0.30,0.40,0.43 in the 
top, middle and bottom rows, respecti vely . The dashed line 
represents P(t) calculated via formula (|ll|). 
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